Simple Maximum-Likelihood Decoding of 
Generalized First-order Reed-Muller Codes 

Kai-Uwe Schmidt and Adolf Finger 



in 
O 

o 

< 



> 



Abstract 

An efficient decoder for the generalized first-order Reed-Muller code RM q (l, m) is essential for the decoding 
of various block-coding schemes for orthogonal frequency-division multiplexing with reduced peak-to-mean power 
ratio. We present an efficient and simple maximum-likelihood decoding algorithm for RM g (l, m). It is shown that 
this algorithm has lower complexity than other previously known maximum-likelihood decoders for RM 9 (1,to). 
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I. Introduction 

IN ORDER to ensure tight control of the peak-to-mean power ratio (PMPR) in orthogonal frequency- 
division multiplexing (OFDM) systems, various block-coding schemes have been proposed [1], [2]. 
These codes are obtained from unions of, say M, cosets of a g-ary generalization of the first-order Reed- 
Muller code RM 9 (l,m) (see Definition [T]). Such codes have the potential to perform error correction and 
ensure a substantially reduced PMPR compared to uncoded transmission. 

Most of the existing decoding algorithms for the above mentioned codes are based on the supercode 
(N) decoding method, described in [3] for the binary case. Such a decoding scheme involves M decodings of 
O ;RMq(l,m). Hence an efficient decoder for RM 9 (l,m) is required. Moreover, if a maximum-likelihood 
^ (ML) decoder is used to decode RM g (l, m), then the supercode decoder performs ML decoding too [3], 

There exist various decoding techniques for the code RM 2 (1, m). The fast Hadamard transform (FHT) 
q can be used to obtain an ML decoder [4]. Suboptimal techniques include majority-logic decoding [4] and 
\ decoding based on the interpretation of the Reed-Muller codes as general concatenated codes [5], [6]. In 
[7] the latter approach was extended to a list-decoding scheme, and it was shown that ML decoding of 
RM 2 (l,m) is possible if the list length is equal to 2. 

In [1] an iterative decoder for RM 2 k(l,m) was proposed that relegates the decoding of RM 2 />(l,m) to 
h decodings of RM 2 (1, m), for which efficient methods exist. A quite similar approach has been reported 
in [8]. Another suboptimal decoder was proposed in [9], where the majority-logic decoding method was 
extended to the nonbinary case. Implicitly, a decoder for RM g (l,m) was obtained in [10], by applying 
the method of ordered statistics to nonbinary codes. In [11] an ML decoder was given for RM 4 (l,m), 
by treating RM 4 (1, m) itself as a union of 2 m cosets of RM 2 (1, m + 1). A g-ary equivalent of the binary 
FHT has been reported in [8] and [12]. Consequently an ML decoder for RM g (l,m) for arbitrary q 
was obtained. In this letter we present a new ML decoder for RM 5 (l,m), which appears to have lower 
complexity than the above mentioned ML decoding schemes. 
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II. Definitions and Preliminaries 

Definition 1: The generalized first-order Reed-Muller code RM 9 (l,m) (m > 1) consists of the Z q - 
valued codewords of length 2 m {u ■ G m (mod q)\u E Z™ +1 }, where we define a generator matrix G m 
recursively by 

Gm-l G m _i\ ... ^ (1 



G »>=^;:; r;:,j » i,h G >-\ i 

Here 0„ and l n are row vectors of length n containing only zeros and ones, respectively. 

The above definition implies the following recursive construction of the codewords of RM ? (l,m) 

u ■ G m (modg) = (u' ■ G m _i|u' • G m _i + u m l 2 m-i) (modg). 

Here (•]•) means concatenation and u = (u ui ■ ■ -Um) and u' = (u Ui ■ ■ -u m _i). This property will be 
used in the next section to derive an efficient ML decoding scheme for RM q (l, m). 

After encoding, the Z^-valued codewords are mapped onto a g-ary phase- shift-keying (PSK) constel- 
lation. This means that for each Z^-valued codeword there exists a corresponding polyphase codeword, 
which is given by 

Here £ = exp(j2n/q) is a primitive qth root of unity, where j 2 = —1. Since there is a one-to-one 
correspondence between Z 9 -valued and polyphase codewords, in the remainder of this letter, we shall 
drop the distinction between them. The context should make clear to which one we refer. 

III. Decoding of RM g (l, m) 

We assume that the codeword £" Gm is transmitted over an additive white Gaussian noise (AWGN) 
channel. Then we receive 

y = C Gm + n, 

where n e C 2 ™ is the complex white Gaussian noise vector. An ML decoder now finds the codeword x 
that is closest to y in the Euclidean sense. We shall call x the ML codeword of y. That is to say 



x = arg mm | \y — x 

i6RM,(l,m) 



|2 



where 

2 m -l 

|2 



\y-x\\ 2 = E 



Vk - x k \ 



k=0 

It is straightforward to verify 

\\ y - x \\ 2 = \ \y\\ 2 + ||a;|| 2 - 2Re{y ■ x H }, 

where () H denotes Hermitian (conjugate transpose), such that y ■ x 11 is the dot product of y and x. Notice 
that \\y\\ 2 and ||a;|| 2 are independent of the decoding result, since ||a;|| 2 = n for each x e RM ? (l,m). 
Hence finding the codeword that is closest to y is equivalent to finding the codeword for which the real 
part of the dot product with the received vector is maximized. To be precise 

x = arg max Re{y ■ x H }. 

xeRM q (l,m) 

We now state our decoding algorithm. Then we show that this algorithm always outputs the ML 
codeword of y. 

Algorithm 2: Soft-decision maximum-likelihood decoder for RM q (l, m). 

1) Input a vector y = (y y 1 ■ ■■y 2 m- 1 ) E C 2 "\ 

2) If m — 1, output the hard decisions x, the estimated information symbols ii, and Re{y ■ x H }. Stop. 



3) For i = 0, 1, • • • ,q — l calculate 



z(i) = (zo(i) zi(i) ■ ■ ■ z^-i-i («)) 

with 

= Vk + Vk+2^C\ 

and use the algorithm to decode z{%) to z(i), to get the estimated information symbols u(i), and 
to calculate p(i) = Re{z(i) ■ z H (i)}. 



4) Determine 



% = argmaxp(i). 



Output the decoded codeword x = (z(i)\z($)£ % ), the estimated information symbols ii = 
and p{i). 

Theorem 3: Algorithm El is a soft-decision maximum-likelihood decoder for RM 9 (l,m). 

Proof: We will show that, if Algorithm E] is an ML decoder for RM g (l, m — 1), then it is also an 
ML decoder for RM g (l,m). The case m = 1 serves as the induction anchor, since it is obvious that 
Algorithm El is an ML decoder for RM 9 (1, 1). Moreover it is clear that, if m = 1, Algorithm El outputs 
Re{y ■ x H } correctly. 

Now let m > 1. We have to show that the algorithm always outputs the ML codeword x and Re{y-x H }, 
which is potentially required for the higher stage of the decoder. From the discussion in SectionHUwe know 
that a polyphase codeword x e RM g (l, m) may be expressed as x = (w\w where w £ RM g (l, m — l) 
and i E Z 9 . By computing 

(w, i) = arg max Rej-y • (w\w C) H \, (1) 
we can obviously find a; = (iu|w;^ 1 ). Let us inspect the real part of the dot product in (HJ) 

2 m-l_ 1 

Re{y(w\wC) H } = J] Re{(y k + y k+2m -i^)w* k } = Re{z(i) ■ (2) 

fc=0 

where ()* denotes complex conjugation and is as calculated in Step 3 of AlgorithmEl By hypothesis 
Algorithm El is an ML decoder for RM 9 (l,m — 1), which is now used to get z(i) according to 

z(i) = arg max Re{z(i) ■ w H }. (3) 

tceRM,(l,m-l) 

Moreover the algorithm provides p(i) = Re{z(i) ■ (z(i)) H }. In order to find w, and thus also x, it remains 
to compute % = argmax i6Z9 p(i), which is done in Step 4 of AlgorithmEl Then we have w = z(i) and i, 
from which x can be obtained. We also have, by and ©, p(i) = Re{z(i) ■ (z(i)) H } = Re{y ■ x H }. ■ 
We close this section with a brief discussion. The code RM 9 (1, m) may be interpreted as a union of q 
cosets of the code {(w\w) \ w e RM 9 (l,m — 1)} with coset representatives {i ■ (0 2 ™-i|l2 m - 1 ) | i G %* q }- 
Then Algorithm El basically performs the steps of the supercode decoding principle, as stated in [3] for 
binary codes. We also mention the relationship of Algorithm El and the list-decoding scheme, proposed in 
[7]. In this reference it was shown that an ML decoder for RM 2 (r, m) can be obtained when the list length 
is equal to 2. The basic idea behind Algorithm El can therefore also be interpreted as a generalization of 
the list-decoding method for nonbinary codes, where the list length is set to q. 



IV. Complexity Analysis and Comparisons 



We next analyze the complexity of Algorithm |2] and make comparisons with existing approaches. In 
practice it is likely that q is a power of 2. So we restrict our analysis to the case q = 2 h . We consider the 
number of complex additions N + {2 h ,m) and the number of complex multiplications N x {2 h ,m) as com- 
plexity measurements. Notice that we regard one real addition as half a complex addition. Multiplications 
with 1, j, — 1, — j are not counted, since they are trivial and can be implemented by sign bit changes and 
swapping of the real and imaginary part. 

We start with m = 1. Then, in Step 2, we need to determine the hard decisions and the dot product of 
two vectors of length 2. The hard decision for the kth symbol is Xk = where 



For h < 2 the above calculation is a trivial task and can be accomplished with simple sign logic. For 
h > 2 we may implement this maximum calculation as follows. First determine the quadrant in which the 
received value lies. This leaves 2 h ~ 2 + 1 candidates, over which the maximum in © has to be calculated. 
For each candidate one complex multiplication is required. Since one of the possible signal points lies on 
the real and another on the imaginary axis (where the multiplication is trivial), we need for each of the 
two hard decisions 2 h ~ 2 — 1 complex multiplications. Hence 



Having Re{y^ Uk } for k = 0,1, one real addition is required to calculate the dot product in Step 2. 
Hence 



Now let m > 1. In Step 3 we require 2 h 2 m ~ 1 complex multiplications to calculate Zk(i) for k = 
0, • • • , 2 m_1 — 1 and i = 0, ■ ■ ■ , 2 h — 1. If h < 2, these are trivial multiplications. For h > 2, 4 out 
of the 2 h complex multiplications for each k are trivial. Moreover, if we have computed Dk£,~ % for i = 
0,1, - ■ ■ , 2 h ~ 1 — 1, by sign bit changes, the values yk£,~ l for i = 2 h ~ 1 , ■ ■ ■ ,2 h — 1 can be easily derived. 
After having the products yk£,~\ there are 2 h 2 m ^ 1 additions required in Step 3. These are real additions 
if h — 1 and complex additions if h > 1. Finally we need to perform 2 h decodings of RM ? (l,m — 1). 
Hence in total we have 



Table HI lists the number of complex operations required to decode RM 9 (l,m) for different m and 
q = 2 h . It also compares the complexity of our algorithm with the complexity of [12, Algorithm 1]. 
Notice that for h = 1 only real additions are required, and its number is equal to half the indicated 
value. It is apparent that Algorithm |2] has lower complexity than [12, Algorithm 1]. To be precise, our 



ii k = arg max Re{y k • £ '}. 



(4) 




N + (2 h ,l) = 1/2. 




and 




TABLE I 

Number of complex operations required to decode RM, (1, m) 
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[12, Algorithm 1] 
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21824 




16 


182528 


449376 
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2 


176 
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4 


4480 





8064 







8 


103680 


87360 


349440 


174720 




16 


2920960 


7190208 


19173888 


14380416 



algorithm requires exactly one half of the multiplications required for [12, Algorithm 1]. For h > 1, the 
ratio of the number of complex additions required for Algorithm |2] and [12, Algorithm 1] is approximately 
{2 h+l + 2 h ~ 1 — l)/2 2h and converges to this value as m approaches infinity. We finally remark that the 
ML decoder for RM 4 (1, m) in [11] requires (m + l)2 2m+1 real additions. Hence it has higher complexity 
than Algorithm El 

V. Conclusion 

We have presented a simple maximum-likelihood decoder for the first-order generalized Reed-Muller 
code RM 9 (l,m) that enables efficient decoding of certain OFDM codes with strictly bounded PMPR. 
It was shown that the presented algorithm has lower complexity than previously proposed maximum- 
likelihood decoding schemes for RM ? (l,m). 
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